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ABSTRACT 



We present results from extensive molecular dynamics simulations of collapse transitions of 
hydrophobic polymers in explicit water focused on understanding effects of lengthscale of the 
hydrophobic surface and of attractive interactions on folding. Hydrophobic polymers display 
parabolic, protein-like, temperature-dependent free energy of unfolding. Folded states of small 
attractive polymers are marginally stable at 300 K, and can be unfolded by heating or cooling. 
Increasing the lengthscale or decreasing the polymer- water attractions stabilizes folded states sig- 
nificantly, the former dominated by the hydration contribution. That hydration contribution can 
be described by the surface tension model, AG = ^(T)AA, where the surface tension, 7, is 
lengthscale dependent and decreases monotonically with temperature. The resulting variation of 
the hydration entropy with polymer lengthscale is consistent with theoretical predictions of Huang 
and Chandler (Proc. Natl. Acad. Sci., 97, 8324-8327, 2000) that explain the blurring of entropy 
convergence observed in protein folding thermodynamics. Analysis of water structure shows that 
the polymer-water hydrophobic interface is soft and weakly dewetted, and is characterized by en- 
hanced interfacial density fluctuations. Formation of this interface, which induces polymer fold- 
ing, is strongly opposed by enthalpy and favored by entropy, similar to the vapor-liquid interface. 
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L INTRODUCTION 

Hydrophobic interactions are one of the major contributors to biological self-assembly in solu- 



ding and aggregation, micelle and membrane formation, and biomolecu- 
5]. Recent work in this area has focused on the lengthscale dependencies 



tion, including protein fo 
lar recognition QjJ, [2, 

of hydrophobic hydration and interactions [4,6, 7, 8,^]. In particular, a recent theory by Lum, 
Chandler, and Weeks (LCW) [6] highlighted the different physical mechanisms of solvation of 
small and large hydrophobic solutes in water. Small solutes are accomodated in water through 



1 1| | - whereas solvation of larger solutes requires forma- 



molecular-scale density fluctuations f | 1 OL 
tion of an interface similar to that between a liquid and a vapor (4j, [6J, LL2D • This change in physics 
is also reflected in thermodynamic (entropy vs enthalpy dominated hydration) [|90 and structural 



(wetting vs dewetting of the solute surface) 
between larger hydrophobic solutes in water 



12, 



1311 aspects of hydration. Similarly, interactions 



MM 11a- 

labdi. 
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1811 are characteristically distinct from 



those between their molecular counterparts H 

The differences between the hydration and interactions of small and large solutes characterize 
many-body effects in hydrophobic phenomena. Effects of similar origin are also at work in as- 
sociation of small hydrophobic solutes into a larger aggregate 112 lL 12211. and are quantified by the 



n-particle potential of mean force (PMF) Q23 



24 
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26J]. For n > 3 however, the dimensionality 



of the system makes calculations of n-particle PMFs computationally prohibitive. 

To this end, ten Wolde and Chandler [27] performed studies of a hydrophobic polymer in water, 
a many-body model system of potential relevance to protein folding. Their calculations, using a 
coarse-grained model of water, showed that a sufficiently long hydrophobic polymer readily col- 
lapses in water driven by a significant free energy difference {^30k B T for their model polymer). 
Their simulations also predicted that the essentially wet extended states of the polymer go through 
dewetted transition states at a smaller radius of gyration before collapsing into an ensemble of 
compact globular states. Although molecular simulation studies of conformational equilibria of 
solvophobic oligomers in explicit water and in other solvents have been performed in the past 



L28, 



32, 



330, a systematic study of how attractive interactions and lengthscales affect 



the thermodynamics of collapse or folding-unfolding transitions has not been reported. 

Motivated by the work of ten Wolde and Chandler Q2VD, here we present results from extensive 
molecular dynamics simulations that address the dependence of folding of hydrophobic polymers 
in water on lengthscales and attractive interactions. We show that compact folded states of smaller 
polymers with attractive polymer-water interactions are only marginally stable and display thermo- 
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dynamic characteristics of both warm and cold denaturation. Increasing the polymer lengthscale 
makes the hydration contribution to folding more favorable, and stabilizes folded states. Reduc- 
ing the polymer- water attractions also stabilizes folded states. We explore the suitability of area 
and volume based models to describe the hydration contribution to folding. Entropy obtained 
from temperature dependence of the hydration contribution to folding varies with the polymer 
lengthscale, explaining the origin of blurring of the 'entropy convergence' in protein folding using 
arguments proposed by Huang and Chandler [Q]. Structural signatures, such as water density and 
its fluctuations in the vicinity of the polymer, indicate that the polymer-water hydrophobic inter- 
face is soft and weakly dewetted and allows enhanced density fluctuations, collectively inducing a 
collapse into compact folded conformations. 
RESULTS AND DISCUSSION 

To observe collapse transitions potentially important to understanding protein folding thermo- 
dynamics, the model polymer needs to be sufficiently long, such that it can adopt compact globular 
as well as extended states. Preliminary studies indicate that 6-, 8-, or 12-mer chains are too short 
to form compact states. In contrast, somewhat longer, 25-mer hydrophobic oligomers, can fold 
into tightly wound helical states. We therefore focus on the collapse transitions of such 25-mer 
hydrophobic chains in water. 

To study the lengthscale dependence, we studied two versions of the 25-mer, a smaller length- 
scale polymer (denoted C25) comprising methane-sized monomers (m) (with <7 mm =0.373 nm and 
m — m bond length of 0. 1 53 nm, and m — m — m bond angle of 1 1 1 .0 deg), and a larger lengthscale 
polymer (denoted CG25) comprising larger (ethane- sized) monomers (M ) (with a MM =0.44 nm, 
and M — M bond length of 0.25 nm). The M — M — M bond angle potential was not employed 
in simulations of the larger lengthscale CG25 polymer, and dihedral potentials were not employed 
for either polymers. We note that we are not interested in detailed studies of specific alkanes in 
water, but in conformational transitions of flexible hydrophobic chains driven by hydrophobic in- 
teractions, and their dependence on lengthscales, attractions, and temperature. The inclusion of 
additional intrapolymer interactions (e.g., a dihedral potential) changes the conformational pref- 
erences but will not affect the qualitative hydration contributions or their temperature dependence 
that are of primary interest here. Further, the larger lengthscale polymer, CG25, can be thought 
of approximately as a coarse-grained equivalent of a 50-mer of methane-sized monomers. The 
M — M — M angle potential was therefore turned off to capture the higher conformational flexi- 
bility of that polymer. 
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To quantify effects of polymer-water interactions on collapse transitions, we employed three 
different interactions for C25 and CG25 polymers: (i) Full LJ, where monomers interact with other 
nonbonded monomers and water via Lennard Jones (LJ) interactions [a mw = (a mm + a ww )/2 and 
£mw = \/€mm£ww\- We used e mm = 0.5856 kJ/mol and e A /M=0.85 kJ/mol, respectively. The 
strength of attractions employed here is similar to that of typical alkane-water attractive interac- 
tions [13411 . (ii) Half LJ, where monomer- monomer e mm was reduced to half its value in (i), and 
monomer-water e mw correspondingly by a factor of \J~2, and, (iii) WCA, where Weeks -Chandler- 



Andersen H35I1 repulsive potential was used to describe monomer-monomer and monomer-water 
interactions. We denote the full LJ and WCA versions of the two polymers by C25-LJ, C25- 
WCA, CG25-LJ, and CG25-WCA, respectively. 

Figure 1 shows the potential of mean force (PMF), W(R g ), for conformational sampling of 
polymers along the radius of gyration reaction coordinate obtained from umbrella sampling sim- 
ulations (see Methods). Several interesting features are apparent. For the C25-WCA polymer, 
tightly wound helical configurations (with R g ~ 0.42nm) are stabilized by ~20 kJ/mol relative to 
extended configurations in the unfolded basin. Although the PMF decreases monotonically as the 
polymer folds, the slope is not uniform - PMF is relatively flat beyond R g of 0.65 nm and rapidly 
decreases below that value. 

For the CG25-WCA polymer comprising larger monomer units, the driving force for folding 
is significantly higher. The folded states are stabilized by almost 100 kJ/mol relative to extended 
states. Indeed, in absence of the umbrella potential, the extended states of this polymer collapse 
rapidly into compact states over a timescale of ~100 ps. Once folded, no significant partial or 
complete unfolding events are observed indicating the strength of hydrophobic driving force for 
the polymer of this size. 

To what extent do attractive interactions change the folding-unfolding of hydrophobic polymers 
in water? Figure 1 also shows PMF profiles for C25-LJ and CG25-LJ polymers interacting with 
full LJ potential having magnitude similar to typical alkane-water attractions. Here, the attrac- 
tive interactions directly affect two opposing enthalpic terms - the monomer-monomer attractions 
favor folding of the polymer into compact states in absence of other effects (Fig. 2c), whereas 
monomer-water attractions favor unfolding (Fig. 2b) to increase the favorable energetic interac- 
tions with water. 

For the smaller, C25-LJ polymer, balance of these interactions leads to two distinct minima 
corresponding to folded and unfolded ensembles separated by a barrier. Overall, the polymer- 
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water attractive interactions dominate over other contributions, leading to a slight stabilization 
of the unfolded states compared to the folded ones. In contrast, for the larger CG25-LJ polymer, 
despite the significant polymer-water attraction contribution favoring the unfolded states (Fig. 2b), 
the qualitative nature of the PMF is similar to that of CG25-WCA polymer - folded states are 
favored significantly over the unfolded states. 

The total PMF, W(R g ), governs the overall conformational behavior of polymers in water, and 
can be separated into three contributions J^ll : 

W(R g ) = W vac (R g ) + (U PW (R g )} + W hyd (R g ), (1) 

where W vac is the PMF profile in vacuum, (Upw) is the ensemble averaged polymer- water in- 
teraction energy (for given R g ), and we define Whyd as me hydration contribution. We note that 
[W(R g ) — W vac (R g )) is traditionally referred to as the 'solvent contribution'. However, its sep- 
aration in Eq. 1 into (Upw(Rg)) and Wh y d{R g ) is similar to that used in perturbation theory 



approaches applied previously to solvation phenomena IIIOL I37I1 . The rationale for writing Eq. Q] 



is that each individual term may be separately predicted using theoretical arguments or simplified 



models [36], as we illustrate below for the Wh y d(R g ) term. For nonpolar solutes in water, the 
solute- water attractive energy is generally smaller in magnitude than the hydration contribution, 
whereas for polar or ionic solutes, it will dominate the solvation free energy. In either case, a 
similar separation may provide a useful predictive framework. 

Figure 2 shows Wh yd (Rg), (Upw{R g )), and W vac (R g ) contributions to W(R g ) for LJ and WCA 
versions of the C25 and CG25 polymers. For both the polymers, the hydration contribution is 
large, approximately 40 kJ/mol and 140 kJ/mol, respectively, and monotonically favors folding 
of polymers into compact states. Despite the significant attractive interactions for LJ polymers, 
Wh y d makes the dominating contribution to the overall PMF. More importantly, the Wh yd {Rg) 
contribution appears to be similar for all C25 or CG25 polymers (full LJ, half LJ, or WCA), 
relatively insensitive to the strength of attraction. This supports the separation of the PMF in 
Eq. 1, and indicates that the Whyd term captures the essential physics of hydrophobic hydration 
that may be described either using cavity formation (for small lengthscales) or interface formation 
process (for larger lengthscales). 

Once the overall W(R g ) profile is available from simulations, the total free energy of unfolding 



a polymer in water, AG U , can be obtained by integration of that PMF 113 811 . Figure 3 shows 



the temperature dependence of AG U , calculated in this manner. Effects of both lengthscale and 
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attractive interactions on the free energy of unfolding are apparent. For the smaller and attractive 
C25-LJ polymer, the AG U profile is parabolic with compact states being only marginally stable 
(by ~1 kJ/mol) for 350 < T < 450 K. Thus, this polymer shows signatures of both warm and 
cold denaturation processes. Increasing the lengthscale increases the hydrophobic driving force 
for collapse. Indeed, folded states of CG25-LJ polymer are stabilized by over 10 kJ/mol at the 
temperature of maximum stability. Reducing the polymer-water attractions or removing them 
completely also increases the driving force for folding. For example, at 300 K, the AG U values 
for C25-WCA and CG25-WCA polymers are about 10 and 30 kJ/mol, respectively. 

The hydration contribution, AG^ yd , to the total free energy of unfolding estimated using 
Whyd(Rg) is shown in Figure 4a. Several features of these curves are noteworthy. The AG^ yd 
is positive at all temperatures, is larger in magnitude for the CG25 polymers, and decreases mono- 
tonically with increasing temperature for all polymers. Because the polymer- water attractions are 
not included in this term, the AG^ yd (T) profile is similar and favors folding of both WCA and 
LJ polymers. Figures 3 and 4 together show that for the purely repulsive WCA polymers, AG^ yd 
determines the temperature dependence of the overall AG U , whereas for the attractive LJ versions 
of the polymer, the additional polymer-water energy contribution (which favors unfolding) makes 
the AG U (T) profile parabolic. 

The hydration contribution quantifies the difference in hydration free energy of folded and un- 
folded state ensembles, AG^ yd (T) = G hyd (T) - G h f yd (T). Thus, to predict AG hyd (T), we can (I) 
explore the applicability of simpler models that predict free energies of hydration of representative 
folded and unfolded states using lengthscale-dependent physics of hydration; or (II) attempt to 
model the entire conformation-dependent potential of mean force W hyd (R g ), and obtain AG^ yd 
by appropriate integration. Below we explore path (I) first. 

The folded globular states of C25 and CG25 have equivalent radius [(3V747r) 1/3 or (A/4tt) 1/2 ] 
of approximately 0.65 and 0.95 nm, respectively; where V is the solvent-excluded volume (SEV), 
and A is the solvent-accessible surface area (SASA) of the polymer. These radii values are in 
the 'crossover' region [|4j, |9J], and we expect the hydration free energy of the folded state to be 



adequately described by the surface tension model, G hyd = jA |36j] . In contrast, the lengthscale 
of unfolded polymers is more difficult to quantify. If we assume the relevant lengthscale to be 
small, equal to that of a monomer, then the free energy of hydration of unfolded polymers will be 
roughly proportional to the solvent excluded volume (4j,|7|], approximated as n x G^ d n (T) (where 
n is the number of monomers, and the subscript mon indicates the value for a monomer), or more 
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quantitatively as, V x c(T), where the volume coefficient c(T) = G^ n /V mon (or alternatively, 

= dGj^on/ dV mon ). 

Figure 4b shows the hydration contribution to the free energy of unfolding predicted using the 
'area- volume' model, AG^ yd (T) = V x c(T) - j(T)A. We obtained V and A values using the 
molecular surface package of Connolly [39]. Values of c(T) (~ 150- 170 kJ/mol/nm 3 ) obtained 
independently from test particle insertions of C25 and CG25 monomers into pure water systems, 
are in excellent agreement with that used by ten Wolde recently [36]. The surface tension, 7(T), 
for folded states can be approximated by that of an equivalent sphere. These data are difficult to 
obtain by test particle insertions [9], but can be estimated using the data in Figure 2 in Huang 
and Chandler [7] (as A/i/A), or by using Tolman corrected macroscopic surface tension, 7(T) = 
7oo(T) (1 - 28(T) [R) (where the Tolman length, 6(T) is taken from |4fl], and 7oo(T) is the surface 
tension of water [J41]), or by using 7(T) = [dWh y d/dA] T . Figure 4b shows the prediction using 
the area- volume model, where 7(T) = [dWh y d/dA] T was used. The prediction not only differs 
numerically from simulation values of AG^ yd , but shows opposite temperature dependence. We 
find that the other two methods of estimating 7(T) give similar results. 

Why does the physically motivated area-volume model not provide a good description of 
AG^ yd obtained from simulations? This model attempts to predict the temperature-dependent 
excess chemical potentials of hydration (or the vacuum-to-water transfer free energies) of folded 
and unfolded states independently - a truly challenging task. We note that only the difference 
AG hyd (T), and not these independent values are currently available from simulations. There- 
fore, we can only make qualitative arguments here regarding the applicability of the area-volume 
model to the present systems. For example, the unfolded states of hydrophobic polymers comprise 
monomers joined together by covalent bonds, as distinct from fully separated and well dispersed 
monomers in solution. Also, there is conformational diversity in the unfolded states that includes 
fully as well as partially extended conformations with different local curvatures. A single small 
monomer-like lengthscale therefore does not appear to describe the hydration free energy of un- 
folded states, and especially, its temperature dependence. In contrast, many other processes do 
involve assembly of small lengthscale moieties into a larger lengthscale aggregate {e.g., in micelle 
formation |42|]) where the area- volume model may provide a good description of hydrophobic 
contributions to the thermodynamics of assembly. 

We also explored the suitability of alternate models (i) AG^ d {T) = c(T)AV u and (ii) 
AG^ yd (T) = j(T)AA u , where AV U and AA U are the volume and areas of unfolding, respec- 
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tively, determined from our simulations. Here, volume or area-based description alone is assumed 
to apply to both folded and unfolded state ensembles. We find that the volume-based calculation 
[model (i)] that uses c(T) obtained from thermodynamics of small solute transfer significantly 
underpredicts the AG^ yd at room temperature, as well as its temperature dependence (not shown). 
In contrast, the area-based description [model (ii)] provides a quantitative prediction of AG^ yd (T) 
as shown in Figure 4a. 7(T) constitutes a critical input to this model, and was obtained using 
7(T) = [dWh y d/dA]T for both polymers, The numerical value of 7 is larger for CG25 polymer 
compared to that for the smaller lengthscale C25 polymer. Temperature dependence, 7(T), for 
both polymers is qualitatively similar to that of vapor liquid interfacial tension measured in exper- 



iments or in simulations of the SPC/E water 114 ill . 

In contrast to the area-volume model, the models (i) or (ii) used above based on volume or 
area alone do not attempt to predict the excess chemical potentials of folded or unfolded states 
independently, but implicitly assume that the PMF Why a varies linearly with volume or area over 
the entire R g range. How good is the assumption of linearity? Figures 4c and d show Wh y d plotted 
as a function of V and A, respectively. All values are referenced to their respective values at R c g at . 
For the CG25-WCA polymer, variation of Wh y d is not linear but roughly sigmoidal with respect to 
V. The slope, c = OWh y d/dV is more than twice that of c mon obtained from small solute studies 
(~ 171 kJ/mol/nm 3 ) and approaches c mon only at the ends. For the CG25-LJ polymer, the behavior 
is similar, if somewhat less non-linear. As the effective lengthscale of the polymer is reduced, the 
slope decreases, and for C25-LJ polymer approaches the c mon (298K) value. Inspection of these 
curves over a range of temperatures shows a similar behavior. For C25 polymers, even though 
W hy d is approximately linear with V at all temperatures, the slope c(T) is different from c mon (T), 
explaining the failure of volume-based AG^ yd (T) = c mon (T)AV u model (i) above. 

Figure 4d shows that Whyd varies roughly linearly with area. The slope, 7 = dWhyd/dA, 
is higher for CG25 polymers and decreases (from 50 dynes/cm to 25 dynes/cm) as the effective 
lengthscale is reduced from CG25-WCA to CG25-LJ to C25-WCA to C25-LJ. The temperature 
dependence of 7 for C25 and CG25 polymers is apparent from Figure 4a (as AA U values are 
constant). 7 decreases monotonically over the temperature range of interest qualitatively similar 
to the vapor-liquid surface tension of water. In absence of detailed simulation data, how one may 
predict a priori the exact numerical value of 7(T) for a polymer of a given lengthscale is not 
entirely clear. Prediction using the Tolman equation (typically used for spherical solutes) provides 
a qualitatively correct picture, but is not quantitative (not shown). Application of LCW theory 
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270 to explore lengthscale and temperature dependence of 7 for non- spherical solutes may 



be one computationally efficient alternative. Nevertheless, the approximate linearity with area and 
the corresponding temperature dependence of the slope explains the sucess of AG^ yd = 7(T) AA U 
model. 

The free energies AG U (T) and AG^ yd (T) obtained here for hydrophobic polymers have im- 
portant connections with the thermodynamics of protein folding. The AG^ yd (T) is positive, fa- 
vors folding, and decreases in magnitude with increasing temperature, a trend similar to that of 
AG U (T) for WCA polymers. The polymer-water attractions in contrast favor unfolding, and more 
so at lower temperatures. Inclusion of attractions therefore reduces the stability of folded states 
and makes the AG U (T) curve roughly parabolic in T, similar to that for proteins. To summarize, 
increasing the lengthscale increases the magnitude of AG U (T) and decreases the temperature, T* 
of maximum AG U , whereas addition of attractive interactions has the opposite effect. 

These observations have two important consequences. First, the increased stability of folded 
polymers with reduction of polymer-water attractive interactions as well as the shift of T* to 
lower values with increasing lengthscale will make it difficult to observe cold denaturation over 
temperature range accessible to experiments, without assistance from other denaturing influences 



B4l 



44fl . Second, because the magnitude and temperature dependence of AG hyd (and therefore, of 



AG U ) is sensitive to lengthscales, the corresponding temperature-dependent entropy contributions 
will also vary with polymer (or protein) size. Sufficient variation in protein lengthscales will 
lead to blurring_of the "entropy convergence" observed in thermal unfolding of globular proteins 



,7 
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Figure 4e shows that for a methane-sized solute, the hydration entropy, AS^ yd = 
— (dAG^ yd /dT), is large and negative at room temperature, characteristic of hydrophobic hy- 
dration at small lengthscales J^ . For C25 and CG25 polymers the AS^ yd is positive over the 
entire range of T, with its magnitude being larger for the CG25 polymer. These observations 
are qualitatively consistent with predictions by Huang and Chandler QVU , who argued that such 
lengthscale dependence of entropy explains the blurring of the 'entropy convergence' observed in 
unfolding of globular proteins. Also, the hydration contribution to the heat capacity of unfolding, 



AC^i = T(dAS^ yd /dT), is positive, similar to that for proteins. 

Is the folding or collapse of hydrophobic polymers observed here induced 



of a vapor-liquid like interface? In their study, ten Wolde and Chandler 112711 found that as a 



by the formation 



hydrophobic polymer folds to more compact states, it goes through dry or dewetted transition states 
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with low vicinal water densities. In Figure 5, we quantify the vicinal water density by calculating 
number, N w , of water molecules in the first hydration shell divided by the polymer SASA. These 
local density values are plotted as a function of the polymer SASA and are normalized to 1 in the 
extended states of polymers. 

The normalization factor, equal to the vicinal density in the extended state, is 14.1 for C25- 
LJ, 13.2 for C25-WCA, 13.3 for CG25-LJ, and 12.8 for CG25-WCA polymers (all in units of 
molecules/nm 2 ). For isolated C25 and CG25 monomers, the local surface densities obtained from 
independent simulations are essentially identical for LJ and WCA versions, equal to approximately 
16 and 15 molecules/nm 2 , respectively. Thus, by this criterion, the extended states of polymers 
themselves are already fairly dewetted compared to isolated monomers, highlighting the differ- 
ences between fully separated monomers and unfolded states of polymers. 

Folding of these polymers leads to further dewetting. The extent of that dewetting is rather 
small, ~3 % for the C25-LJ polymer. When the attractive interactions are removed (e.g., for 
the C25-WCA case), the local density decreases, goes through a minimum (~10% reduction) 
before reaching the folded states, for which the reduction is about 6%. This depletion is consistent 
with a small increase in the lengthscale upon folding of the polymer. These effects are enhanced 
for the larger CG25 polymers, which show density depletion even for attractive polymer-water 
interactions (Figure 5b). 

The fluctuations of water density in the hydration shell are also noteworthy. Figures 5c 
and 5d show normalized water number density fluctuations in the first hydration shell, a 2 Nw x 
SASA/ (N w ) 2 , which quantify approximately the compressibility of the first hydration shell. (We 
note that a 2 N also shows similar features). It is clear that for the small attractive C25-LJ polymer, 
the hydration shell compressibility changes little as it folds. In contrast, the larger and repulsive 
CG25-WCA polymer shows increased absolute fluctuations as well as a peak in the magnitude of 
fluctuations as the polymer folds to compact states. Interestingly, although the nature of W(R g ) 
and the dewetting is qualitatively similar for CG25-WCA and CG25-LJ polymers, an explicit peak 
in fluctuations is not observed during folding of the CG25-LJ polymer. This is expected as the 
polymer-water attractions will dampen the longitudinal density fluctuations. In addition, solute- 
water attractions increase the crossover lengthscale, or alternately, reduce the effective lengthscale 
of the solute. In this sense, CG25-LJ lengthscale is smaller than that of CG25-WCA polymer. 

Collectively, for the polymers of the nanometer lengthscale studied here, we find that the un- 
folded states are dewetted relative to separated monomers, and the folding process is accompanied 
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by further dewetting of the polymer surface. More importantly, the increased density fluctuations 
for repulsive polymers indicate that the polymer-water interface is soft, qualitatively similar to 
water- hydrophobic liquid-liquid [50] or liquid- vapor interfaces. In addition, the positive hydration 
entropy indicates that the unfavorable free energy is dominated by the enthalpic contribution sim- 
ilar to that for interface formation [4, ^], supporting the above assertion. Extension of the above 
results would predict increased dewetting, higher density fluctuations, and larger driving forces 
for the folding of longer polymers. 
CONCLUSIONS 

We have presented results from extensive molecular dynamics simulations of folding-unfolding 
of hydrophobic polymers in water with focus on understanding effects of polymer lengthscale and 
polymer- water attractive interactions on thermodynamic and structural aspects of polymer folding. 
The folding of model polymers provides a route to characterizing the many-body hydrophobic ef- 
fects in the context of realistic biological self-assembly process, such as protein folding. We 
studied conformational transitions of 25-mers comprising monomers of methane-like and larger 
particles. To study the effects of polymer- water attractions, we employed attractive LJ as well 
as purely repulsive WCA interactions between monomers and water. We found that the smaller 
attractive polymer displays a parabolic protein-like unfolding free energy, AG U (T), profile as a 
function of temperature. The folded states of this polymer are only marginally stable at room 
temperature and can be unfolded by heating or cooling the solution. Increasing the lengthscale 
or reducing the polymer-water attractions not only increases the stability of folded states, but also 
shifts the temperature of maximum stability to lower values. Also, the dependence of hydra- 
tion entropy on polymer lengthscale is consistent with predictions of Huang and Chandler [7|] on 
lengthscale dependence of hydrophobic hydration that provide a basis for understanding of the 
blurring of entropy convergence observed in protein unfolding thermodynamics. 

The hydration contribution to folding of hydrophobic 25-mers studied here is large and strongly 
favors compact folded states. We explored the applicability of physically motivated models based 
on area and/or volume descriptions of hydration of polymer conformations to predict the hy- 
dration contribution, AG^ yd . Our analysis indicates that predictions using AG^ yd = ^{T)AA U 
are in excellent agreement with simulation data. The surface tension 7 is larger for CG25 poly- 
mers and decreases monotonically with temperature for all polymers. Structural details indicate 
that the polymer-water interface, especially in the vicinity of the folded state ensemble, is soft, 
weakly dewetted, and is characterized by enhanced density fluctuations. Correspondingly, AG^ yd 



12 



is enthalpy-dominated consistent with interface-driven folding process. 

The polymer models presented here or the ones used by ten Wolde and Chandler, present 
promising model systems of potential importance in understanding the role of hydration in thermo- 
dynamics of protein folding and aggregation [51]. The ability to change intrapolymer and polymer- 
water interactions, and polymer size allows for systematic studies of the free energy landscape in 
the extended parameter space. Further, investigation of the folding-unfolding of these polymers 
over a range of solution conditions, such as high pressures, in presence of cosolutes and solvents 
113 ll B2I1 . will provide important insights into protein stability over the broader thermodynamic 



space. 

METHODS 

Intra polymer interactions: The intra polymer potentials for bond lengths in C25 and CG25 
polymers and bond angles in C25 polymers were harmonic in nature with following parameters: 
For m — m bond length in C25 polymers, V& = 0.5k b (r — r ) 2 , where k b = 334720.0 kJ/mol/nm 2 
and r = 0.153 nm, and for m — m — m bond angle, Vg = 0.5ko(9 — #o) 2 , where kg = 462.0 
kJ/mol/deg 2 and 9 = 1 1 1.0 deg. Similarly, for M — M bond length in CG25 polymers, we used 
k b = 60702.0 kJ/mol/nm 2 and r = 0.25. 

Free energy calculations: We used the umbrella sampling technique |52|] to characterize the 
conformational free energy of hydrophobic polymers in water as a function of polymer radius of 
gyration, R g . A harmonic potential, \k u (R g — R g o) 2 , was applied to efficiently sample polymer 
conformations in a window near R g0 . For C25 polymers, we used 13 equally spaced windows 
spanning a range of R g between the most compact (~ 0.37 nm) and the most extended (~ 0.92 
nm) states, whereas for CG25 polymers, we used 20 equally spaced windows with R g ranging from 
~ 0.45 nm to ~ 1.50 nm. Simulations of C25 polymers included 1700 water molecules (SPC/E 
model Il53ll ). For CG25 polymers, depending on R g0 value, the number of water molecules ranged 
from 4000-12000. Additional simulations in which 12000 water molecules were used uniformly 
in all windows confirmed that the results presented here are insensitive to system size. 

The values of k u ranged between 4000 and 12000 kJ/mol/nm 2 , based on the local gradient 
of the free energy, dW(R g ) /dR g , estimated from preliminary simulations. Data from simulation 
trajectories - 26 ns for the C25 polymer and 40 ns for the CG25 polymer (2 ns/window) - were 



combined using the WHAM formalism [54, 



55 



5611 to generate a potential of mean force profile 



along the R g coordinate for each polymer. Similar calculations were performed in vacuum to 
obtain the PMF in the absence of solvent. 



13 



5811 , with suitable 



Simulation details: MD simulations were performed using GROMACS 11571 . 
modifications to implement umbrella sampling as well as the WCA interaction scheme. Periodic 
boundary conditions were applied and the particle mesh Ewald method [590 was used to calculate 
the electrostatic interactions with a grid spacing of 0. 1 nm. NPT simulations were performed at 1 
atm and 300 K maintained using Berendsen algorithm [60]. The SETTLE algorithm JotI was used 
to constrain OH and HH distances in water with a geometric tolerance of 0.0001 A. Simulations 
were also carried out over a range of temperatures to calculate temperature dependence of the free 
energy of folding along the saturation curve of water [47]. 
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FIG. 1: The lengthscale- and attraction-dependent hydrophobic collapse. The PMF, W(R g ), at 298 K, for 
the smaller C25 and larger CG25 polymers interacting with water via attractive (LJ) and repulsive (WCA) 
interactions. Conformations with R g < R c g ut are defined as compact states: R c g ut = 0.6 nm for C25, and 
0.73 nm for CG25 polymers, respectively, as shown by arrows. PMF for C25 and CG25 polymers are 
zeroed at R g values of 0.83 nm and 1.4 nm, respectively. To make the comparison clear, the horizontal axis 
is R g - R™ in ; with R™™ = 0.385 nm for C25 and =0.474 nm for CG25 polymers, respectively. 

FIG. 2: The hydration and other contributions to PMF for polymer folding at 298 K. (a) Wh y d{Rg) con- 
tribution for C25 and CG25 polymers interacting via full LJ (filled circles), half LJ (inverted triangles), 
and repulsive WCA (open circles) interactions. Panel (b) shows the polymer-water attractive contribution, 
(U pw (R g )}. Lines are guide to the eye. Panel (c) shows the intra polymer W vac (R g ) contribution for LJ and 
WCA versions of C25 and CG25 polymers. 

FIG. 3: Free energy of polymer unfolding. The free energy AG U obtained by integrating W(R g ) curve as 

IJut cxp(—W(Rg)/k B T)dR g 

exp(— AGu/ksT) = ^ cui is shown by circles. The curves are obtained by fitting 

S R Ln eM-W(R g )/k B T)dR g 

the data to form AG U {T) = AH u (T ref ) - TAS u {T ref ) + AC P _ U [(T - T ref ) - T\n(T/T ref )] shown by 
solid lines, where T re f = 298 K was used. 
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FIG. 4: Hydration contribution to the free energy and entropy of polymer unfolding, (a) AG^ yd (T) es- 
timated from simulations (symbols) by subtracting appropriately integrated polymer-water and intrapoly- 
mer energy and entropy contributions: LJ polymers (filled circles), WCA polymers (open circles), with 
red and blue indicating C25 and CG25 polymers. Intrapolymer contributions were obtained from vac- 
uum runs of the polymers at the same temperature. Lines in panel (a) are predictions of the area model, 
^Qhyd _ j(t)AA u , where AA U for C25 and CG25 are 0.73 and 1.66 nm 2 , respectively, obtained by 
taking differences of ensemble average area of folded and unfolded states, and j(T) = [dWh y d/dA] T . (b) 
Comparison of AG^ yd predicted using the area-volume model (lines) for C25 and CG25 models with sim- 
ulation data (symbols). Panels (c) and (d) show the relative conformational hydration free energy, AW/,^, 
plotted as Wh y d{R g ) — Wh y d(Rg Ut ), as a function of solvent excluded volume (SEV) and SASA, respec- 
tively. SEV and SASA are plotted relative to their average values for Rg Ut conformers, indicated by AV 
and AA. Dashed lines in panels (c) and (d) with respective slopes c(298 K) and 7(298 K) (values indi- 
cated) bracket the observed behavior for different lengthscales. Panel (e) shows the hydration contribution 
to entropy of unfolding. Entropy of hydration for WCA-methane obtained using test particle insertion 



calculations [37] is also shown 



FIG. 5: Water density fluctuations and weak dewetting of polymers, (a) and (b) The local water den- 
sity in the vicinity of C25 and CG25 polymers, respectively, calculated as average number of hydration 
waters divided by average SASA, (N w ) /SASA. These values are normalized to 1 in the wet extended 
states of the polymers, (c) and (d) A measure of hydration shell compressibility or fluctuations obtained 
as 0% / {N w ) 2 SASA, where crjy is the variance of water number fluctuations in the hydration shell of 
the polymer. A water molecule is considered to be in the hydration shell if it is less than a certain cut off 
distance (0.6 nm for C25 and 0.7 nm for CG25) from any monomer of the polymer. 
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Figure 1,Athawaleet al. 
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Figure 2, Athawale et al. 
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Figure 3, Athawale et a I. 
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Figure 4, Athawak 
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Figure 5, Athawale 
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